1 Introduction 



Both the matrix eigenvalue problem and the inverse eigenvalue problem are 
the subjects of intensive research [1,2]. This paper addresses both problems. 
A means of producing matrices with any specified eigenvalues and whose 
eigenvectors are known is presented. Both the elements of the matrix and the 
eigenvectors of the matrix are given in analytic form; this should facilitate the 
use of such matrices for analysis in other areas of linear algebra. A reversal of 
the matrix-generation procedure naturally suggests a matrix-diagonalization 
algorithm, and we derive the equation which needs to be solved to realize 
this. The findings we are going to present grew out of research into quantum 
foundations [3], but are here presented from a purely mathematical point of 
view. 

The organization of this paper is as follows. In the next section, we 
present the basic result that makes possible the production of matrices with 
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any desired eigenvalues. In Section 3, we outline the diagonalization routine 
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that arises from it. In Section 4, we give the functions which for the case 
ijT) • N = 5 enable the practical realization of both these projects. The results of 

preliminary application of this theory are reported in Section 5. We conclude 
^ with brief comments in Section 6. 
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2 Matrix Generation 

Suppose that the iV 2 functions <j>{Bi\ Cj) possess the orthonormality property 
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1 £<l>(B l ;C i )<t>*{B l ;C j ) = 6 ij , (1) 
i=i 

where B and C are parameters that take the values Bi,B2,...,Bn and C\, 
C2,...,Cjv respectively. The quantities (j)(Bi]Cj) are functions of some vari- 
ables, collectively denoted by x. If with the aid of the N arbitrary numbers 
Ai, A2, ...Ajv we form the N x N matrix [M] using the rule 
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M l3 = <f>*(B i; C n )<P(B f C n )X n (2) 

71=1 

we find that the eigenvectors of this matrix are 
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with the respective eigenvalues Aj. Indeed if we write 

[v} = MU 

we find that the kth row of [M] [£J is 
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= \4>*(B k -Ci). (5) 
Thus the eigenvalue equation 

Mfo]=Aifo] (6) 
is satisfied. Property Eqn. (0) means that the eigenvectors are orthonormal: 



5ij. 



(7) 

These results give us a means of generating matrices with specified eigen- 
values and known eigenvectors. We only need appropriate functions 4>(Bi, Cj). 
An example of such functions will be given in Section 4. Given such func- 
tions, and a set of eigenvalues {A,;}, we are able generate an infinite number 
of different matrices by varying the argument x of these functions. Each such 
matrix is generated together with its own orthonormal set of eigenvectors. 

On the other hand, given a value of x, we can generate an infinite number 
of matrices by using different sets of eigenvalues. In this case, all these 
matrices share the same set of orthonormal eigenvectors. 
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3 A Diagonalization Algorithm 



Since each matrix is generated by means of a specific value of the argument 
x, we should be able, at least in principle, to find a way of deducing this 
value of x for a given matrix. A procedure for doing this would constitute a 
diagonalization routine. In this section, we derive one such method. 

Suppose we are given a matrix [M] formed by the prescription Eqn. (J2J) . 
Its eigenvectors are given by Eqn. ©. If we only knew the value of x that 
was used to form the matrix, then we could fix the eigenvectors and recover 
the eigenvalues. In order to develop a way of deducing the value of x, let us 
start, for the sake of convenience, by setting 



so that 



^C i ;B j ) = (f>*(B j ;C i 



(8) 



*(B 2 ; C, 



(9) 



V £(Ci,B N ) J 



Now, from 



and the value of x that was used to produce the matrix, we can recover the 
eigenvalue Aj from the rth row by using 



(A, 



Y,i Mri£(Ci, Bi 



(11) 



£(Cj, B r ) 

The sth row would give the same value. Hence, comparing the values 
from the rth and sth rows, we can write 

{Xi)r ~ {Xi)s = 0. (12) 

This must be true for all possible combinations of r and s : 
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In other words: 
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(14) 



But if the value of x used is not correct, the equation is not satisfied. 
Thus, this equation is a means of deducing the correct value of x. A form of 
this equation better adapted to numerical work is 
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aC t , B s ) £ M rl Z(Q, Bi) - £(Ci, B r ) £ M sl C(C l7 Bi) 
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= 0. 



(15) 

We have summed over all the eigenvalues by means of the index i because 
we must ensure that the equation is satisfied for every eigenvalue. 

Once the correct value of x (denoted by x ) is found, the eigenvalues are 
calculated using 



(K)r — 



£,(Ci, B r ) 



(16) 
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and the eigenvectors are 



(17) 
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As is wise in the solution of such equations, it is necessary to establish 
that the value of x obtained is not spurious by putting it back into the 
equation. In this case, this means using it to calculate the eigenvalue using 
every line of the eigenvector. Thus, for each eigenvalues, there should be iV 
values, which should agree in order for this value of x to be acceptable. Of 
course, in the case where one of the elements of the eigenvectors vanishes, it 
cannot be used to compute the eigenvalue, and there will be fewer than iV 
values to compare. 



4 Probability Amplitudes 

In order to utilise the prescription for generating matrices with specified 
eigenvalues and to realise the algorithm for diagonalizing them, we need the 
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functions 4>(Bf, Cj). Since the origins of the method presented in this paper 
lie in quantum theory, we shall call these functions probability amplitudes. 
A source for such functions is spin theory in quantum mechanics[3]. The 
treatment of a spin s system yields functions suitable for generating matrices 
of order N — 2s + 1. As an example, we take the spin-2 system. In the 
direction of the unit vector 

c = (sin 9 cost/?, sin 9 sin ip, cos 9) (18) 

its operator is [4] 



[c-S] = 
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(19) 

The eigenvalues of this matrix are 2, 1, 0, —1 and —2, with the respective 
normalized eigenvectors 
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(20) 
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be 
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and 
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(24) 



/ 



We have used the index m to label the eigenvalues of the spin matrix 
because that is the convention. Also, we need to keep the distinction between 
the eigenvalues of this matrix and the eigenvalues {Aj} belonging to the 
matrix we wish to generate once we have obtained the probability amplitudes. 

For the direction 



b = (sin 9' cos ip', sin 9' sin <p' , cos 9') 



(25) 



the eigenvectors of the operator are identical, except that (9, ip) is replaced 
by (9',ip'). We now associate the parameter B with the direction b, and 
the parameter C with the direction c. The notation is such that the val- 
ues £>i, I?2, B 5 correspond respectively to the eigenvalues 2, 1, — 2 . The 
same holds for C . 

Some familiarity with spin theory in quantum mechanics is perhaps nec- 
essary to fully appreciate these functions. Each element in each eigenvector 
is a probability amplitude, and its modulus square is a probability for a 
certain measurement [3] . However, from a mathematical point of view, all 
that is needed is the fact that these elements possess the property Eqn. (JTJ). 
The scalar products of the eigenvectors corresponding to b with those corre- 
sponding to c give the most generalized forms of these functions. Thus the 
required quantities are [5] 



^B l ;C j ) = [ X ^[x i S}. 

Evidently, in this case, x = (9, <p, 9' , p'). 

To illustrate the notation, we give a few of these functions: 



(26) 
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[xQ=2]^[Xm=2] 

cos 4 - cos 4 -e* 2 ^ 
2 2 

+ sin 0' sin 6> cos 2 - cos 2 ^e^"^ + | sin 2 0' sin 2 6> 
2 2 8 

+ sin 0' sin 9 sin 2 - sin 2 -e^'*') 
2 2 

. 4 . 4 6* _j2((fl-(p') 

+ sin - sin — e ^ ^ ; , 



(27) 



and 



sin 9 cos 4 - cos 2 -e* 2 ^') 
2 2 

i/o ■ 2 ^ 2® \ ■ qI 2 ® 2 ® i(v-v') 

+ 3 sin cos - ) sin 9 cos — cos - e KV v ' 
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3 

— sin 2 9' sin cos 6* 
4 

-(3 cos 2 - - sin 2 -) sin 9' sin 2 - sin 2 B -e*W> 
v 2 2 J 2 2 

-sin0sin 4 ^sin 2 ^e- 2 ^') (28) 



0(B 5 ;C 5 ) = [xE- 2 ] f [x£-2] 



4 ^ 4 ^ i2(w-w') 

= cos - cos — e w ^ ; 
2 2 

+ sin 0' sin 9 sin 2 - sin 2 -e^-v 7 ) + - s i n 2 0' sin 2 9 
2 2 8 

+ sin 9' sin cos 2 ^ cos 2 9 -e^~^ 

+ cos 4 ^cos 4 ^e^ 2 ^'). (29) 
There are 22 other functions. We remind ourselves that since 

£(C i ,B j ) = <l>*(B j ;C i ), (30) 
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then 

m,B J ) = ix$nx®]- (si) 

These are the functions which enter into Eqn. (j!5|) . 

In general, the probability amplitudes satisfy the inter-dependence law 



[Bi, Cj) 
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where d is a third direction to which corresponds the new parameter D. 
This is a particular form of a fundamental quantum property of three sets of 
probability amplitudes belonging to one quantum system [6] , 
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Y,x{A i ;B l ) ( j>{B l ;C i 



(33) 
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In the context of quantum theory, the quantity \4>{Bi\ Cj)\ 2 is the proba- 
bility that if the spin projection in the direction b is rriih, a measurement of 
it in the direction c yields the result rrijh[3]. Here Ms Planck's constant. 

The parameters B and C, as observed, are defined by the angles (#' , <p') 
and (9, ip), and the indices of each parameter correspond to the eigenvalues. 
When (9', ip') = (9, ip) in Eqn. (0 , A = C and Eqn. {TJ) is satisfied. Eqn. © 
is a special form of Eqn. (J33|) . Within the context of quantum measurement 
theory, all this has a ready interpretation [6]. However, we stress that it is 
not necessary to understand this in order to use these functions to generate 
matrices customized as to eigenvalues and type. 

With these probability amplitudes, the eigenvector for the eigenvalue Aj 

is 



/ 



V 



'(B 2 ;a) 

'(-B 4 ; Ci) 
'(B s] Ci) 



I [xL b i 2 ] t 

(b) 



[Xm=l 

w itr 



(34) 



] [xQ] 

V [xE-2] f [x@] J 

Of course, since the matrix Eqn. (JT§j) is a particular case of the form of 
Eqn. (j2J), its eigenvectors are given by Eqn. evidently, they correspond 
to the angles ^' = and <p' = 0. In other words, for this case, the unit vector 
b is in the z direction. 
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The association between the eigenvalues of the matrix we are generating 
and their eigenvectors is through Eqn. (J2J), where the value of C and the 
eigenvalue have the same index. 

5 Results 

A Fortran 77 program written to generate the matrices confirms that dif- 
ferent families of matrices are obtained depending on the values of 9, tp, 9' 
and p' ,a.nd on the character of the eigenvalues {Aj}. The following observa- 
tions are made, and can be predicted from the structure of the probability 
amplitudes and of Eqn. (J2J). 

If all the angles are zero or if 9 = 9' and tp = tp', the resulting matrix is 
diagonal. 

If the eigenvalues are real, the matrix is Hermitian, but its elements 
depend on what values of the angles are used. 

Whatever the eigenvalues, the matrix is symmetric if tp — tp'. 
If tp — tp', the eigenvectors are real. 

If all the eigenvalues are pure imaginary, and the arguments are arbitrary, 
the matrix is anti-Hermitian. 

If all the eigenvalues are pure imaginary, and tp = tp', the matrix is imag- 
inary but symmetric. 

If the eigenvalues are arbitrary and the values of the arguments also 
arbitrary, the matrix is general. 

The kinds of families that can be generated with different combinations 
of the eigenvectors and the arguments have not been fully investigated, but 
it seems probable that it does not require much ingenuity to generate such 
special forms as tridiagonal matrices, etc. Clearly, with this method, much of 
the inverse-eigenvalue problem is solved. Further investigation should show 
how to chose the values of the angles so as to obtain the desired matrices 
with specific values of certain elements. 

As far as the diagonalization algorithm is concerned, it has only been 
partially validated at this time. The main complication here is that the non- 
linear equation (|15j) is in four variables. Therefore, solving it is somewhat 
complicated. The following was however achieved. A matrix was generated 
by means of Eqn. (J2J). In order to diagonalize it, three of the arguments were 
held at their actual values, and the remaining one was determined with the 
use of Eqn. (|15p. Which one of the variables to treat as unknown could be 
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decided at pleasure. The value of this variable was determined from Eqn. 
(JTSJ) by means of the bisection method. That this approach was successful 
indicates that once a fast and reliable method of solving a non-linear equation 
in four variables is employed, the algorithm will prove to be of some utility. 

6 Conclusion and Discussion 

In this paper, we have presented a prescription for forming matrices in such 
a way that their eigenvalues and eigenvectors are known. The method is 
very general indeed, and simply by varying the arguments, different kinds 
of matrices can be obtained. Always, the normalised eigenvectors are simul- 
taneously given. The eigenvectors can be chosen to be real if desired, and 
the any combination of eigenvalues can be used. By essentially reversing the 
matrix-generation procedure, we have proposed an algorithm for diagonaliz- 
ing matrices formed this way. 

A good amount of work is still needed in order to fully understand and 
utilize the methods presented. For example, it is necessary to classify more 
properly and completely according to the values of the arguments 8, tp, 9' and 
<// the kinds of matrices that can be generated. Such information would be 
particularly helpful in the solution of the non-linear equation. If a matrix is 
such that one or more of the arguments must have certain values, this makes 
so much easier the job of solving the equation, since the number of variables 
is reduced. 

The matrices dealt with here are of order 5. In order to deal with matrices 
of higher order, it is necessary to have the functions {<f\ for those cases. A 
source of these functions will always be spin theory, but treatment of other N- 
dimensional quantum systems should produce other forms of the functions. 
Each such set of functions probably produces matrices of different characters. 
As such, it expands the range of uses to which we can put these matrices. 
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